Giant direct magnetoelectric effect in strained multiferroic heterostructures 
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The direct magnetoelectric (ME) effect mediated by lattice strains induced in a ferroelectric film 
by a ferromagnetic substrate is evaluated using first-principles-based calculations. To that end, the 
strain sensitivity of ferroelectric polarization and the film permittivity are calculated as a function 
of the in-plane biaxial strain for Pb(Zro.52Tio. 43)03 films under various depolarizing fields. It is 
found that the ME voltage coefficient varies nonmonotonically with this strain and may reach giant 
values exceeding 100 Vcm^^ Oe~^ over a strain range that can be controlled through the electrical 
boundary conditions. 



Multiferroic materials exhibiting magnetoelec- 
tric (ME) effects represent the subject of cut- 
ting-edge research since they are promising for 
applications in the next generation of microelec- 
tronic and nanoelectronic devices^. As the room- 
temperature ME effects in single-phase materials 
are weak, the theoretical and experimental stud- 
ies are currently focused on composite multifer- 
roics. Here the interface-related and proximity ef- 
fects may strongly increase the ME response, as 
found, in particular, for the strain-mediated ME 
eflects in ferroelectric-ferromagnetic hybridal. 

The design of advanced multifunctional ma- 
terials may be greatly facilitated by theoreti- 
cal modeling. In particular, ah initio and first- 
principles-based calculations already proved to be 
able to predict physical properties of single crys- 
tals and thin films directly from microscopic phe- 
nomena occurring at the atomic scale^. Applica- 
tion of first-principles calculations to ferroelectric- 
ferromagnetic heterostructures and composites, 
however, is a challenging task as the behavior of 
such systems is much more complex than that of 
single-phase materials. Recently, several studies 
were performed on multiferroic heterostructures 
with the 2-2 connectivity^"—. However, the pre- 
dicted ME effects are either weak or extremely 
localized at the interfaces, which limit their pos- 
sible applications in microelectronic devices. In 
contrast, the strain-mediated ME effect spans over 
the whole volume of the active phase in a properly 
designed heterostructure. Accordingly, the macro- 
scopic ME response becomes strong as well, which 
is critical for the development of ME sensors of 
magnetic fields and many other devices^. 

Theoretical predictions of the strain-mediated 
direct ME effect generally require the determina- 
tion of three ingredients: (i) the deformations in- 
duced in the ferromagnetic component by the ap- 



plied magnetic field; (ii) the degree of strain trans- 
mission through the interface between two ferroic 
constituents of a hybrid material system; and (iii) 
the response of polarization and dielectric suscep- 
tibility of the ferroelectric phase to changes in lat- 
tice strains. In the case of epitaxial heterostruc- 
tures, the mechanical coupling between two phases 
is usually very strong so that the strain transmis- 
sion may be perfect. The deformation response of a 
ferromagnet to the measuring ac magnetic field (5H 
becomes large only in the presence of a bias mag- 
netic field H2rJ^. The dependence of this response 
on H is nonmonotonic and cannot be calculated 
theoretically with sufficient accuracy^^. Accord- 
ingly, the first ingredient should be evaluated using 
the available experimental data on the field depen- 
dence of magnetostrictive deformations. On the 
contrary, the properties of a single- crystalline fer- 
roelectric phase subjected to strains can be calcu- 
lated from first principles with a high accuracyi^. 

In this Letter, we calculate the strain-mediated 
direct ME effect for a multiferroic hybrid in the 
form of a thin ferroelectric film sandwiched be- 
tween two electrodes and coupled to a thick 
ferromagnetic substrate. This geometry maxi- 
mizes the ME coefficients since the film does not 
suppress the magnetic-field-induced deformations 
of the substrate. As a representative example 
of a ferroelectric material, we have chosen the 
Pb(Zro.52Tio.48)03 (PZT) disordered solid solution 
with the morphotropic composition. The strain 
sensitivity of the out-of-plane polarization and the 
film permittivity are calculated for (OOl)-oriented 
PZT films by first-principles-based methods. Since 
epitaxial films may be strained to a different ex- 
tent during the deposition, these characteristics are 
calculated as a function of initial biaxial in-plane 
strain r]„i typically induced in perovskite films de- 
posited on (OOl)-oriented cubic substrates. The 



bottom electrode is assumed to be fully strained 
by the substrate or thick buffer layer due to co- 
herent lattice matching at the interface. The pres- 
ence of electrodes is taken into account in the cal- 
culations through the j3 parameter that describes 
partial screening of the depolarizing fieldi^. The 
film is modeled by a 12x12x12 supercell with peri- 
odic boundary conditions along the [100] and [010] 
pseudocubic directions (chosen as the xi and X2 
axes, respectively). The nonperiodic [001] direc- 
tion (x3 axis) thus allows for a finite thickness of 
the film. 

The effective Hamiltonian used in our calcula- 
tions is described in Refs.— . Importantly, it takes 
into account not only the ferroelectric local modes 
{ui) (i — 1,2,3) and strains (f^y), but also rota- 
tions of the oxygen octahedra (w^) as well as the 
fluctuations of m^, Ty^j, and uji. This feature, to- 
gether with the account of the depolarizing-field 
and finite-thickness-related effects, represent the 
major differences from the previous phenomenolog- 
ical calculations^^. As a result, our study provides 
a significantly improved description of the ME ef- 
fect, notably bringing into light the crucial role of 
the depolarizing field. 

Mechanical boundary conditions imposed on a 
thin epitaxial film involve three fixed components 
of the strain tensor, which for the cube-on-cube 
epitaxy become 7711 = 7722 — rim and 7712 = 0. The 
biaxial misfit strain rjm equals rjm — afUm/aref^^, 
where afum is the actual in-plane lattice parameter 
of the film and are/=4.065A is the lattice constant 
of the prototypic cubic phase of bulk PZT at room 
temperatursi^i. It is worth noting that afUm coin- 
cides with the lattice parameter of the substrate in 
the case of coherent epitaxy, but may differ from 
it considerably in hybrids with misfit dislocations. 

The total energy of the supercell was used in 
Monte-Carlo (MC) simulations with 4.10'' MC 
sweeps. To compute the strain-mediated ME ef- 
fect, the code was modified to run calculation as a 
function of strain at constant temperature (300 K). 

Figure [1] displays the predicted effect of biax- 
ial strain on the polarization components Pi cal- 
culated using the standard relationship between 
Pi and Ui and the out-of-plane dielectric con- 
stant 633 under partially screened depolarizing field 
{(3—0.98). With the strain rjm varying from com- 
pressive to tensile, the equilibrium ferroelectric 
phase changes from the tetragonal c one (space 
group PAmm) to the monoclinic r phase (Cm) 
and then to the orthorhombic aa phase {Cmm2). 
The permittivity €33 strongly increases at the r-aa 
phase transition because P3 goes to zero here. The 
account of rotations of oxygen octahedra changes 
the film polarization behavior dramatically, de- 
spite the fact that the mean angle ((w)) is found 
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FIG. 1: Polarization components (a), out-of-plane di- 
electric constant (b) and average magnitude of the fluc- 
tuations of the local oxygen octahedra tilting in each 
unit cell (c) as a function of biaxial misfit strain for a 
48A-thick PZT film at 300 K under partially screened 
depolarizing field (/3 = 0.98). The influence of oxy- 
gen octahedra rotations (labelled AFD) is illustrated 
on panels (a) and (b). 



to be zero, in agreement with refji^. Indeed, the 
decrease of P3 with increasing compressive strain 
is linked to the average magnitude of the fluctua- 
tions of the local oxygen octahe dra ti lting in each 
unit cell (see the evolution of \/(uj^ with 77^ on 
FiglTJc)), as when rotations are frozen P3 increases 
slighty (see FigjIJa)). In addition, these rotations 
change the position of the r-aa phase transition 
and associated dielectric peak and affect the max- 
imum value of permittivity 

The film strain states leading to the formation 
of various phases shown in FiglT] can be achieved 
through an appropriate buffer layer deposited on a 
ferromagnetic substrate. When the magnetic field 
H is turned on, the film strains change due to 
substrate magnetostrictive deformations. In the 
symmetric case, where H is orthogonal to the film 
plane and the substrate has a four-fold symme- 
try about the X3 axis, these changes reduce to a 
variation of the biaxial strain 77,„. Hence the ME 
effect can be evaluated via the strain sensitivity 
dP^/drjm of the out-of-plane polarization P3 with 
respect to 77™ . 

The magnetic field parallel to the film surfaces, 
however, breaks the initial isotropy of in-plane 
strains. To quantify the ME effect appearing in 
this situation, one has to calculate the sensitivi- 
ties of P3 with respect to uniaxial strains 7711 and 



7^22- Therefore, we computed the evolutfon of the 
local mode U3 as a function of uniaxial strain rjaa 
superimposed on a biaxial strain yy^. Represen- 
tative results obtained for the three ferroelectric 
phases are shown in Figl2j It is seen that rjn and 
7^22 have the same effect on M3, which is consistent 
with the symmetry of the phases. In the c and 
aa phases, the calculated dependences are linear. 
In contrast, the evolution of the local mode occur- 
ring in the r phase follows a power law. However, 
on the scale of small strains induced by the mea- 
suring magnetic field ^H the linear approximation 
proposed in Ref,-'^ holds. It should be noted that 
the dependences shown in Figl^] also enable us to 
allow for the infiuence of the bias magnetic field 
H (applied along the [100] or [010] direction). The 
strain changes induced by H, however, usually may 
be neglected in comparison with initial strains be- 
cause relevant magnetostrictive deformations are 
well below 10-3^-J^. 
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FIG. 3: Sensitivities of the out-of-plane polarization to 
7711 (squares), 7)22 (diamonds), and r/m (with AFD full 
circles, without AFD empty ones) plotted as a function 
of the initial biaxial strain -qm- The lines are cubic 
spline fits to the data points. 
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FIG. 2; Out-of-plane local mode U3 as a function of 
uniaxial strain t^h or 7722 (closed and open symbols, 
respectively) in the (a) c, (b) r, and (c) aa phases. 
The lines represent the best fits with a linear ((a) and 
(c)) and power-law (b) functions. 



The strain sensitivities dPs/drjaa calculated 
from the local mode U3 are plotted in Figl3] as a 
function of the initial biaxial strain. As expected 
from the symmetries of predicted phases, within 
the accuracy of calculations dP^/drju = dPz/drj22- 
At the same time, these sensitivities are signifi- 
cantly smaller than the sensitivity dP-i/dr]m of po- 
larization with respect to the biaxial strain itself 
(with and without AFD), which is also plotted in 
FigEl This result is consistent with the prediction 
of the thermodynamic theory^ giving dP^/drjm = 
2 dP^/drjii owing to the superposition principle. 
The positions of the maxima, however, are shifted 
away from the r-aa transition, in contrast to ther- 
modynamic theory. This is caused by inhomoge- 
neous strains arising at the local scale within the 
film because strains are calculated in each cell of 



The calculated strain sensitivities and dielectric 
constant enable us to evaluate the polarization 
{ctp^j—dPs/dHj) and voltage {aE'ij=dE^/dHj) 
ME coefficients characterizing the discussed multi- 
ferroic hybrid. Since the figure of merit of a mag- 
netic field sensor is the output voltage it delivers, 
we focus on the ME voltage coefficients govern- 
ing the electric field i?3 induced in the ferroelec- 
tric parallel-plate capacitor. If the magnetic field 
is parallel to the [100] or [010] crystallographic axis 
of the film and does not induce shear strains here, 

we have as3,— ^ (|^d™i + f^d™,), where 
djaa ^'I'S the substrate piezomagnetic coefficients 
at a given bias magnetic field^'^. Taking the max- 
imum longitudinal coefficient (i5^]^=d^2— 5 10~^ 
mA"^ as measured for the FeBSiC alloy (see Fig. 
3 in refi^) and estimating the transverse coeffi- 
cient to be d722=rf"ii=-2.510-®m k'^ from the 
condition of volume conservatioi>i^, we calculated 
the ME voltage coefficients a^sj as a function of 
initial strain rj„i . Figure S] shows that within the 
accuracy of calculations aE'ii—aEZ2 and further 
demonstrates a nonmonotonic strain dependence 
of these transverse ME coefficients. Remarkably, 
they reach a giant value of ~150 Vcm^^ Oe~^ in 
the middle of the stability range of the r phase, 
which significantly corrects the prediction of the 
thermodynamic theor y ^^1^^ . Moreover, aE^i and 
Q^B32 change sign near the c-r phase transition, i.e. 
where dP^/drjaa does. Accordingly, the ME re- 
sponse vanishes at a certain misfit strain so that 
even the presence of out-of-plane polarization does 
not guarantee the existence of a non-zero ME ef- 
fect. 

The longitudinal ME coefficient aEss can be 
evaluated from the relation aE33=—j--g^dlYi, 



where ^3"! is the effective piezomagnetic coefficient 
characterizing local magnetostrictive deformations 
beneath the film/substrate interface. Since the 
substrate magnetization here is orthogonal to the 
interface, the demagnetizing field is not negligible, 
in contrast to a plate-like substrate subjected to 
an in-plane magnetic field, which was considered 
above. As a result, the relevant local field may be 
much smaller than the applied field, which makes 
'^311 << ^^122 fo'" ^ plate-like substrate. How- 
ever, the effective coefficient ^3^1 can be increased 
strongly by employing a rod-like substrate^ -. Since 
aE33 is not very sensitive to ^3^, we plotted this 
coefficient in Fig|3] assuming that dg^^ = ^"22 . Al- 
though this plot represents the upper bound of 
ctE33, we see that the longitudinal ME coefficient 
may reach giant values about 200 Vcm^^Oe^^ 
even if d^'^i is 20% smaller than d^22 ■ 



SrRuOa electrodes have /3 exceeding 0.99 at the 
considered thickness t=4.8nm. 

Figure [S] shows the longitudinal ME voltage co- 
efficient calculated under different screening condi- 
tions. It can be seen that the highest ME response 
is achieved at zero depolarizing field (/3=1.00). At 
the same time, the maximum value of the ME co- 
efficient remains very high even at /3 = 0.95. Re- 
markably, with decreasing screening ability of elec- 
trodes, the peak of ME response shifts from posi- 
tive to negative misfit strains. This is in agreement 
with the effect of changing /3 on the temperature- 
misfit strain phase diagram of BaTiO^^fl. This fea- 
ture provides a new route for the development of 
ultra-sensitive ME sensors since it shows that the 
ME response can be enhanced not only by tuning 
the misfit strains"'^^ but also by selecting an appro- 
priate electrode material. 
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FIG. 4: Magnetoelectric voltage coefficient of a multi- 
ferroic hybrid comprising a 48A-thick PZT film cou- 
pled to a ferromagnetic substrate subjected to a 
magnetic field parallel to the [001] (bullets), [100] 
(squares), or [010] (diamonds) crystallographic axes 
of the film. The lines are cubic spline fits to the data 
points. 

Finally, we studied how the ME performance 
of the device depends on the screening ability of 
electrodes. This effect can be probed by chang- 
ing the parameter /? that governs the amount of 
screened depolarizing fieldi^. Although formally /3 
may range from unity to zero, for a device with 
metallic electrodes only values close to unity are 
relevant. Moreover, since the calculations showed 
that at /3 < 0.95 the monoclinic r phase cannot 
be stabilized at room temperature, we restricted 
our analysis to 0.95 < (3 < 1.00. It should be 
noted that /? can be estimated from the relation 
/3 = 1 — eo/(eo + Cji), where Ci is the capacitance 
density associated with two film-electrode inter- 
faces, eo is the permittivity fo the vacuum, and 
t is the film thickness. Using the values of Ci de- 
termined by first-principles calculations^^, we find 
that typical ferroelectric capacitors with Ft and 



FIG. 5: Longitudinal ME voltage coefficient calculated 
under different electrical boundary conditions varying 
from the perfect screening of the depolarizing field 
(/3=1.00) to an imperfect screening with 13—0.95. The 
lines are cubic spline fits to the data points. 

In summary, we performed the first cal- 
culation of the strain-mediated ME effect by 
a first-principles-based technique. It was 

found that a multiferroic hybrid comprising a 
Pb(Zro.52Tio. 43)03 thin film coupled to a ferro- 
magnetic substrate may display giant longitudinal 
and transverse ME voltage coefficients exceeding 
100 Vcm~^ Oe~^. Remarkably, this giant ME re- 
sponse exists in a wide strain range Srjm^O.Ol, and 
the ME peak can be positioned exactly at the ini- 
tial biaxial strain r]m existing in the film by an 
appropriate choice of electrodes. 
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